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We compute the Casimir energy between an unusual pair of parallel plates at finite temperature, 
namely, a perfectely conducting plate (e — > oo) and an infinitely permeable one (fj, — > oo) by applying 
the generalized zeta function method. We also compute the Casimir pressure and discuss the high 
and the low temperature limits. 

PACS numbers: 35.20.-i. 



I. INTRODUCTION 

Since Casimir's paper [|l| on the attraction between two 
parallel perfectly conducting plates due to the vacuum 
fluctuations of the electromagnetic field was published, a 
considerable amount of effort, which varies from the in- 
vestigation of new geometries and theories to the applica- 
tion of the Casimir effect to alternative technologies, has 
been put into the study of this importante subject. (For 
a review see, for example, Mostepanenko and Trunov ||] 
or Plunien et at) . Recently, the experimental observation 
of this effect was greatly improved by the experiment 
due to Lamoreaux From the theoretical viewpoint 
Casimir's approach to this problem essentially consisted 
in computing the interaction energy between the plates 
as the regularized difference between the zero point en- 
ergies with and without boundary conditions dictated by 
the physical situation at hand, for instance, the perfectly 
conducting character of the plates. The great novelty! 
indeed of Casimir's 1948 paper was not the fact that two 
neutral object were atracted to each other, which was fa- 
miliar to those studying dispersive van der Waals forces, 
but the simplicity of the method of calculating this at- 
traction in the framework of qu antum field theory. 

Casimir's definition of the vacuum energy requires a 
regularization recipe for its implementation. Many reg- 
ularization techniques are available nowadays and, de- 
pending on the specific physical situation at hand, one 
of them may be more suitable than the others. Par- 
ticularly, methods of computing effective actions are in 
general very powerful to give physical meaning to the di- 
vergent quantities we must deal with. Here we will be 
concerned with one of these methods, namely, the so- 
called generalized zeta function method [Q. We will ap- 
ply it to the case of a pair of parallel infinite plates one of 
which is perfectly conducting (e — > oo), while the other is 
infinitely permeable (/x — > oo). The setup will be consid- 
ered to be in thermal equilibrium with a heat reservoir 
at finite temperature T. This problem at zero tempera- 
ture was analyzed by Boyer || two decades ago in the 
framework of random electrodynamics, a kind of classi- 
cal electrodynamics! whi ch includes the zero-point elec- 
tromagnetic radiation. Boyer was able to show that for 



this unusual pair of plates the Casimir energy is positive 
which results in a repulsive force per unit area between 
the plates. Here we shall study this problem within the 
framework of finite temperature QFT in order to take 
into account the thermal effects in equilibrium. We will 
employ units such that Boltzmann constant, the speed of 
light and Planck's constant are set equal to unity. 

II. EVALUATION OF THE FREE ENERGY 

Since we will be dealing with a system in thermal equi- 
librium, the imaginary time formalism will be convenient. 
In order to apply the generalized zeta function method, 
let us introduce the partition function Z for a bosonic 
theory |: 

Z = N J [D(j>] exp dr J d 3 xcj , (1) 

where C is the Lagrangian density for the theory under 
consideration, TV is a constant and 'periodic' means that 
the functional integral is to be evaluated over field con- 
figurations satisfying: 

<)>(x,y,z,Q) =<j)(x,y,z,0), (2) 

where f3 = T~ 1 , the reciprocal of the temperature, is 
the periodic length in the Euclidean time axis. The 
Helmholtz free energy F(/3) is related to the par- 
tition function Z{j3) through the relation F(f3) = 
— \ogZ(P). Other than the periodic conditions given 
by (0) , we must also consider boundary conditions which 
are determined by the geometry and the nature of the 
physical system under study. An example is the config- 
uration mentioned above. Choosing Cartesian axes such 
that the axis OZ is perpendicular to both plates with 
the perfectly conducting plate at z — and the infinitely 
permeable one at z = d, the boundary conditions on the 
vacuum oscillations of the electromagnetic field are the 
following: the tangential components of the electric field 
as well as the normal component of magnetic field must 
vanish at z = 0, while the tangential components of the 
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magnetic field must vanish at z = d. For the plate ge- 
ometry tha! t we are considering, the electromagnetic 
field can be mimicked by a scalar massless field <f>. The 
boundary conditions stated above can be translated into: 

At n\ n d(j)(T,x,y,z = d) 

<p(r,x,y,z = 0) = 0; = 0, 3 

oz 

where r is the Euclidean time. The insertion at the end of 
the calculation of a factor 2 will take into account the two 
possible transverse polarizations of the electromagnetic 
field. Thus we write \ogZ{0) as: 



log Z = 



-- ) logdet {-d E \T d ) 



(4) 



where d E = d 2 /dr 2 + V 2 , and the symbol Td stands 
for the set of functions which satisfy conditions (j|) 
and (Q). The generalized zeta function method ba- 
sically consists of the following three steps: (i) first, 
we compute the eigenvalues of the operator — d E sub- 
ject to the appropriate boundary conditions and write 
£(s; — d E ) = Tr(— d E )~ s ; (ii) second, we perform an ana- 
lytical continuation of £(s; —d E ) to a meromorphic func- 
tion on the whole complex s-plane; (iii) finally, we com- 
pute det(-d E \T d ) = cxp (- aC( %°J" aE) ) ■ Combining 
equation (^) with the definition of free energy we obtain: 



F(f3) = -[3 



_ 1 d ((s = 0;-d E ) 
ds 



(5) 



The eigenvalues of — d E whose eigenfunctions satisfy 
and (0) are: 
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where k Xl k v £ 5R, n € {0,1,2,3,...} and m € 
{0, ±1, ±2, ...}. The generalized zeta function then reads: 



«s,-8 E ) = L> EE/I 



m— — oo n—0 * 



k 2 x + k 2 y + (2n+l) 2 — + 



4-K 2 m 2 
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(7) 



where L 2 is the area of the plates. After rearranging 
terms in the summations, changing to polar coordinates 
and integrating the angular part out, we can rewrite this 
last equation as: 




(8) 



where k 2 = k 2 + k 2 and the prime on the summation 
symbol serve to remind us that the integer n assumes 



odd values only. Using the following representation for 
the Euler beta function, c.f. formula 3.251.2 jjj: 



f <^+«T-'= fit- 1 -''- f) 



u+2f-2 



(9) 



where 



B[x,y) = 



r(x)T(y) 



(10) 



T(x + y) 

which holds for 3? (v + < 1 and 3?/i > 0, we obtain: 
L 2 r(s - 1) 
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(11) 



In order to connect the simple sum on the r.h.s. of the 
above equation to the Riemann zeta function £r we write: 
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On the other hand, the double sum can be expressed in 
terms of Epstein functions which for any positive integer 
N and 9? z large enough are defined by MM : 
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where oi, ...ajv and M 2 > and writing: 
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we can write: 
C(«, ~d E ) = 



l 2 r(a-i) 
4^ r(s)7r 2 ~ 2s 

2-2s 
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The Epstein functions can be analytically continued to 
a meromorphic function in the complex plane, (see for 
example ||). For N — 2 and M 2 = the analytic con- 
tinuation is given by: 
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Here K v [z) is a Macdonald's function]^] defined on the 
complex z-plane cut along the negative real axis, [— oo, 0] . 
Performing the appropriate substitutions for z, cii and a-i 
and taking advantage of the useful fact that the derivative 
of the function G(s)/T(s) at s = is simply G(0) for a 
well-behaved G(s) we obtain: 
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It follows that the Helmholtz free energy per unit area 
for this peculiar arrangement is given by: 



L 2 ~ 8 X 720d 3 flf ^ In J 
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where we have already accounted for the two possible 
polarization states. The first term in (|l^) represents the 
regularized repulsive Casimir energy at zero temperature 
found by Boyer 0|. Notice that this term is —7/8 times 
the result obtained for the Casimir effect with Dirichlet 
boundary conditions at T — 0. The second term in ( |l8| ) is 
the contribution to the free energy due to thermal effects 
and we can recast it into a more manageable form as is 
shown next. 

The Macdonald's functions K v {z) of half- integral order 
are given by (c.f. formula 8.468 in M): 
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Hence, defining the dimensionless variable £ by £ := 
d/nft = Td/ir and making use of (|l9|), we can recast 
(Eq) into the form: 



*We use the terminology employed by N. N. Lebedev in Spe- 
cial Functions and Their Applications, Dover Publications, 
New York, 1972. The function K v (z) is also known as mod- 
ified Bessel function of the third kind and Bessel function of 
imaginary argument. 



F{p) 7 
L 2 8 



720d 3 irfi 



;/(£), 



(20) 



where /(£) is a dimensionless function defined by the 
double sum: 
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The sum over n can be readly evaluated and after some 
manipulations we end up with: 
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Equation (|22j) summarizes all thermodynamical informa- 
tion concerning the bosonic excitations confined between 
the plates. From ( p2[ ) we can easily obtain the low tem- 
perature regime of the free energy. It suffices to set 

coth (jg) ~ 1, sinh (^j^J ~ CX P (^f) an< ^ keep tne f erm 
corresponding to n = 1: 
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This yields the low temperature limit 
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The very high temperature limit is obtained by setting 

coth f |f) ~ an< ^ s ^ n ^ (st) ~ i| an< ^ evaluating the 
sum. The result is: 
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This will lead to the Stefan-Boltzmann term correspond- 
ing to a slice of vaccum of volume L 2 d. More accurate 
results at high temperature demand that we transform 
the slowly convergent sum over m in (|T]) into a more 
rapidly convergent one. This can be accomplished with 
the help of Poisson summation formula as we shall see 
next. 



III. THE PRESSURE 

Let us go back to the scaled free energy /(£) defined by 
the double sum ([H]). Each term in (|2l|) can be written 
in the form 



' \ m 3 m 2 



-nm/c 



°° poo 

j/3 2 E / dujcj]n(l- 

^, i J 71 K 



(26) 



3 



where a, b and c are constant satisfying the condition 
a = be, and k = 1/ /3c. The scaled free energy can be 
recast into the form 
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where ki = 7r/2d and K2 = 7r/rf. The first term ( |27| ) cor- 
responds to the thermal correction for two infinite par- 
allel perfectly conducting (e — > 00) plates separated by 
a distance 2d. The second one corresponds to the same 
setup but with the plates separated by a distance d. In 
this way we can see that the setup we are considering here 
is the difference between the two setups described above. 
If we apply this reasoning to the zero temperature term, 
see (P0|), we reproduce the factor 7/8 with the correct 
algebraic sign. The net pressure on the plates is given 
by minus the derivative of the free energy per unit area 
with respect to the distance d between the plates and 
like the free energy it splits into the zero temperature 
contribution and the thermal corrections, that is 
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The thermal contribution reads: 
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Now we are ready to make use of one of the several ver- 
sions of Poisson summation formula fl0|| . The particular 
version suitable for our purposes reads: 
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If we use Q30| ) in (|29j) and add the result to the zero 
temperature contribution we obtain 
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This result holds for all temperatures. Notice that the 
zero temperature pressure is apparently missing in our 
final result. This happens because upon the applica- 
tion of Poisson summation formula we obtain, besides the 



Stefan-Boltzmann term and the two sums, a term with a 
negative sign which exactly cancels out the repulsive zero 
temperature contribution. A similar cancellation occurs 
also in the case of the Casimir effect for confined massless 
fermions at finite temperature JHJ and in the high tem- 
perature limit of the standard electromagnetic Casimir 
effect as shown, for instance, in Plunien et al ||. Nev- 
ertheless, it is a straightforward matter to show that if 
we take the zero temperature limit of ([HJ) we recover the 
zero temperature term. 

The high temperature limit is also easily obtained 
from (31). Approximating conveniently the hyperbolic 
cotangents in the sums and evaluating the second partial 
derivatives and keeping the leading correction terms only 
we obtain: 
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A simple integration of (|3l| ) yields another possible rep- 
resentation for the Helmholtz free energy of this setup: 



F((3) 



1 d 1 



7T 2 d 



45/3 4 327r 3 /? 3 <9££ 



coth (47r 2 m£) 



Tfl 



1 d 1 



8tt 3 /3 3 di £ 



E 



coth (27r 2 m£) 



(33) 



where the integration constant is determined by demand- 
ing that in the very high temperature limit the only sur- 
viving term in ( |33"|) must be the Stefan-Boltzmann term. 
Notice that we can also determine this integration con- 
stant analyzing the zero temperature limit of (|33|). The 
high temperature limit of (p3|) is given by 
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Apart from the all-important signs and numerical factors, 
these results compare to those obtained in this limit for 
the attractive case, see for example Plunien et al [E[. 



IV. CONCLUSIONS 

In this paper we have shown how neatly the generalized 
zeta function rcgularization method applies to the repul- 
sive electromagnetic Casimir effect at finite temperature 
for the simple geometry of infinite parallel plates with 
certain special physical properties. As a follow up of the 
application of this method we have obtained expressions 
for the Helmholtz free energy and the force per unit area 
acting on any one of the two plates which comprise this 
peculiar system. We also obtained the low and high tem- 
perature limits of those two quantities. It almost goes 
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without saying that the standard Casimir effect at zero 
and finite temperature can also be treated in the same 
way. 

It is also worth noticing that our high temperature 
limit result is compatible with ideas of dimensional re- 
duction that occurs for T — > oo. Remark that in this 
limit, if we omit the Stefan-Boltzmann term, the domi- 
nant term in the pressure is proportional to l/d 3 , a result 
reminiscent of 2 + 1 dimensions, in contrast to the l/d 4 
term, typical of 3 + 1 dimensions. 

The interesting similarities and differences of the 
Casimir effect associated with a massless bosonic field, 
which arise when we compare the consequences of impos- 
ing Dirichlet boundary conditions in one case and mixed 
ones in the other, indicate that a similar investigation in 
other theories, such as the massive scalar field at zero as 
well as at finite temperature might be rewarding. This 
investigation is being carried out and results will be pub- 
lished elsewhere. 
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